Time correlations in ID quantum impurity problems. 



F. Lesage, H. Saleur* and S. Skorik 
Department of physics, University of Southern California, Los-Angeles, CA 90089-0484- 

(November 28 1995.) 



in 

On 

o 

Q 

(N 



> 
oo 

o 

ON 



o 



We develop in this letter an analytical approach using 
form-factors to compute time dependent correlations in inte- 
grable quantum impurity problems. As an example, we ob- 
tain for the first time the frequency dependent conductivity 
G{uj) for the tunneling between edges in the f = 1/3 frac- 
tional quantum Hall effect, and the spectrum S{ij) of the 
spin-spin correlation in the the anisotropic Kondo model and 
equivalently in the double well system of dissipative quantum 
mechanics, both at vanishing temperature. 

PACS numbers: 72.10.Fk, 73.40. Hm, 75.10.Fk. 



Several problems of great physical interest can be re- 
formulated as impurity problems in one dimensional Lut- 
tinger liquids. These include the anisotropic Kondo prob- 
lem, the double and multi-well problem in dissipative 
quantum mechanics, the tunneling between edges in the 
fractional quantum Hall effect. Although these problems 
are integrable, exact results have for a long time mostly 
concerned static, thermodynamic quantities. A recent 
approach using massless scattering has also allowed the 
computation of DC transport properties and noise, in 
and out of equilibrium. Unfortunately, AC properties 
had so far remained inaccessible, while containing most 
of the open physical questions. In this letter, we de- 
velop an approach based on "massless form-factors" to 
compute the AC properties. We illustrate it on two sim- 
ple cases, both at vanishing temperature: the frequency 
dependent conductivity G{uj) for the tunneling between 
edges in the v = 1/3 fractional quantum Hall effect, and 
the spectrum S{lo) of the spin-spin correlation in the dou- 
ble well problem of dissipative quantum mechanics and 
in the anisotropic Kondo model. 

Edge excitations in the fractional quantum Hall ef- 
fect with filling fraction ly = l/{2p -\- 1) are thought to 
be described by a chiral Luttinger liquid 0. This al- 
lows a number of interesting theoretical predictions to 
be made, in particular concerning the effect of impurities 
For a single impurity (as can be obtained experi- 
mentally through a constriction), one expects, based on 
perturbation theory and scaling arguments, the (dimen- 
sionless) DC conductance to behave as c T^(i/'^~i) at low 
temperature and v-i-c' T'^i'^~^^ at high temperature. The 
problem can be described by the hamiltonian : 

1 

H=— dx[Il^ + {d^<t>f]R + [H^ + {dA?]L 

+ A cos[<^L(0)-(^fl(0)] (1) 



where 0i//j(a;±i) are left and right moving field moving 
on different edges. This model can be simplified through 
the change of variables : 



<l>''{x,t) = —[(j)L{x,t) + (j)R{-X,t)] 

'/'°(a;,t) = -^[(j}Lix,t) - <j)R{-x,t)]. 



(2) 



The even field factorises and after folding the real line we 
get a hamiltonian on the half line with interaction at the 
boundary which can be rearranged in the form : 



H 



1 



IGm/ 



dx [tf + id^(t>f] -f A cos[(/)(0)/2]. (3) 



This hamiltonian is integrable, and the problem is thus 
"solvable" in principle. The DC conductance Q and 
the DC shot noise [D, have been computed in the 
general case of a non vanishing voltage. These results 
followed from combining the standard thermodynamic 
Bethe ansatz with a Boltzmann equation. 

The problem of a two state system in a dissipative bath 
can also be mapped onto a boundary problem . When 
the dissipation is ohmic, it can described by a single spin 
interacting with a bath of electrons. The corresponding 
hamiltonian is the anisotropic Kondo model, which is also 
known to be integrable: 
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where a± are Pauli matrices. The models (||) and (Q) 
have the same "bulk" part: a free boson (we denote the 
coupling i> or a by 5 in the general discussion to fol- 
low), but they differ significantly by the boundary inter- 
action. They both belong to the same large class of in- 
tegrable boundary field theories, and can be approached 
by a unified formalism 0. The general strategy is to 
describe the bulk part by a basis of states that scat- 
ter in a very simple way on the boundary. To find the 
right basis, one can think of more general hamiltonians 
that look like (^ and (^) but with an additional term of 
the form A J_^dx cos[(l){x)]. The bulk part is then the 
well known sine-Gordon model, while it can be shown 
that the boundary interaction does not spoil integrabil- 
ity The right basis for the sine-Gordon model is well 
known, it consists of (massive) solitons/antisolitons, and 
for g < 1/3, [l/g — 2] breathers. These have factorized 
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scattering, and they also scatter simply, without particle 
production, at the boundary. One can now take the limit 
A ^ to obtain a description of the free boson using 
massless solitons/antisolitons and breathers [|o|. This 
is of course a more complicated basis than the standard 
plane waves, but in the latter, the effect of the bound- 
ary interaction is essentially intractable, while here it is 
transparent. Massless particles have dispersion relation 
E = p (resp. E — —p) for right (R) (resp. left (L) ) mov- 
ing particles. Since the theory is massless, we can set the 
arbitrary energy scale equal to unity, and parametrize the 
energies by a rapidity: E = for solitons/antisolitons, 
E = 2sin{mTge^ /2{1 — g)) for the n*'' breather. 

In what follows, we work in Euclidian space and choose 
X to be the imaginary time. This is a "modular trans- 
formed point of view" , where the boundary interaction 
does not appear in the hamiltonian any more, but is 
encoded in a boundary state \B >, so correlators can 
be represented as < OlOOji? >. The boundary state 
has a simple expression [^| in terms of solitons/anti- 
solitons and breathers creation operators, Z*{(3), with 
e = ±,l,...,[l/.g-2] : 

\B >= J2 J K'''''{I3b ~ f3i) . . . K'"''Hf3B - Pn) 

n,€' s 

z;i(/3i) . . . zi^{(3^)z:,j,{(3,) . . . z;, ^(/3„)|0 >, (5) 

where the integrals run on [—00,00] and are ordered 
/?!<...</?„. The boundary interaction is completely 
encoded in the matrix if"'', which derives from the re- 
flection matrix, solution of the boundary Yang-Baxter 
equation: 



(6) 



The strength of the boundary interaction is encoded in 
an energy scale Tb = e^^ ■ The latter is related in a 
non-universal way to the microscopic coupling, Tb oc 
Ai/(i-9). 

To compute a correlation function in this approach, 
one needs to know the matrix elements of the operators 
in the massless particle basis, the so-called form factors 
^t| . Away from g = 1/2, there are infinitely such matrix 
elements to compute, since the theory is truly interac- 
tive. For instance, the U{1) current acting on the vac- 
uum can create arbitrary numbers of solitons / antisolitons 
pairs and breathers. There are thus two difficulties: to 
compute the matrix elements, and to sum their contri- 
butions. In this letter, we will show how to compute the 
form-factors by taking appropriate limits of the massive 
sine-Gordon form-factors. We will also show that expan- 
sions in multiparticle processes converge extremely fast, 
so in practice only a few terms are necessary to obtain 
excellent precision (10^^) on the quantities of interest, 
and this all the way from the short distance to the large 
distance fixed point. 



We start with the frequency dependent conductance 
in the problem of tunneling between edge states. Using 
Kubo formula and folding the system as explained 
above one finds G = f + AG(ti'), where : 

^G{wm) = p2 / dxdx' / dy e'^-'y 

[<d,cj){x,y)d,,(bix',0)> +cc], (7) 

ujm being a Matsubara frequency, z = x + iy. We have 
used here the fact that in the current-current correlation 
the term < dz(j)dz'4> > and its conjugate are insensitive 
to the impurity p^ , contributing | to the conductance. 
This is because dz(j) (resp. dz(j)) act only on R (resp. L) 
particles, while the effect of the boundary is to mix L and 
R-particles. Considering a particular term in the expan- 
sion of \B > with n L- and n R-particles, the only process 
with a non vanishing amplitude in < 0\dz(j}dz'(j)\B > is 
to have dz'4' annihilate all the n L-particles and then 
dz4> all the n remaining R-particles. We thus need the 
form-factors : 



fiPl, Pn)ei,...,<i„ — 



liPn)\0> 



(8) 



Current form-factors have been computed by Smirnov 

Pfor the massive sine-Gordon model. Expressions for 
can then be obtained by taking a massless limit, ie 
by sending the physical mass to zero and the rapidities 
to infinity to keep excitations of finite energy p3[ | . 

As examples, when g = 1/2, there is only a pair of 
solitons/antisolitons in the spectrum. The only non zero 
form factor is /±=p(/3i, /32) = Tc e'^^^'^^^''^^ , with c a 
known normalisation. For g — 1/3, a breather appears 
in the spectrum and the first two form factors are given 
by: 



/i(/3) = cie^ 



for the 1 breather and : 

/±,t(/?1>/?2) = TC2 



C(/3l - 132) g(/3i+/32)/2 

cosh(/3i - P2) 



(9) 



(10) 



for the two solitons form factors. Here C(/3) is a known 
function and ci,C2 known constants, whose expressions 
we will give elsewhere 10 . 

From (||), the form factors expansion results into the 
general expression 

< dz(l){x,y)dz'(l){x',y') >= 

/ dE g{E)exp[E{x + x')-iE{y-y% (11) 
Jo 

and from (^ : 



AG(w) = — Img(-iu;) 
4a; 



(12) 
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Using the previous formulaes and the expression for 
the form factors and the boundary scattering matrices 
IBl, we find : 
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Giuj) = ^[l- iTB/u;)ti,ii-\Lu/TB)] 



(13) 



for 5 = 1/2 in agreement with previous results ||T^. For 
g — 1/3, there is an infinity of form factors contribut- 
ing to the correlation. However, when computing the 
conductivity, we find a very rapid convergence with the 
number of rapidities; two rapidities are sufficient to get a 
1% accuracy. Let us stress that this convergence is inde- 
pendent of the strength of the impurity, and the results 
are valid for the whole flow from small to large distances. 
The quantity G{uj) is plotted in figure 1 as a function of 
Tb/uj. 

Some general features of G{uj) can easily be deduced 
from this approach. The reflection matrices of soli- 
tons/antisolitons expand as a double power series in 
exp f3 and exp(i — l)/3, the reflection matrix of breathers 
as a power series in exp/3. This leads to a double power 
series in (w/Tb)~^+^/3 and (w/Tg)^ at small frequencies, 
(Tb/w)^"^^ and {Tb/oj)^ at large frequencies. Therefore, 
as first argued in ||l^, at low frequencies, G{lu) goes as 
for g < 1/2 and f^j. ^ > 1/2. 

The same method can be applied to compute the spin 
spin correlation G{t) — ^ < [cr(t), (t(0)] > or its fourier 
transform conventionally denoted x"(a;) in the two state 
problem A difficulty arises at first sight because the 
massless scattering description of the anisotropic Kondo 
problem involves only the massless sine-Gordon particles 
and no spins (physically, this is because this description is 
based on the large distance limit of the theory where the 
spin is screened). However, using the fact that spin-flips, 
which are induced by a± , are coupled with insertions of 
vertex operators at the boundary, one can relate C{t) to 
the current current correlator, and get the expression : 
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Im[g{~iij,TB)-g{-tuj,0)], (14) 



where Q is the Fourier transform of the L-R current cor- 
relator deflned in ( [III) . The only difference with the con- 
ductance problem is in the boundary interaction. The 
reflection matrix for the solitons arc now given by ||l0|| : 



Rl = tanh 



and for the breathers we find 



tanh(f 



4(1/9-1) 
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The computations are then done along the same lines 
as before. For g = 1/2, the free or Toulouse point, the 
previously known result : 
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is recovered. For a whole domain of (7 G [.25, .6], the 
form factor expansion gives again very precise results for 
x"(w) with only the first two terms, for all strengths of 
the impurity coupling . Let us give as an example the 



explicit expressions for g 
breather is : 



i . The contribution from the 
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with c = —0.141 and the contribution for the two solitons 
form factors is given by : 



-Tie 



' dpM^zM^^^,, (19) 
00 cosh2(/3-log(l-e^)) 



RtiP + l0g(;^))i?+(l0g[(l - e'');^]) - 1 



'Tb 



Tb' 



with c' — —0.0451 and the function C(/3) given in pi[ . 
In figure 2 we show the function S{uj) = x"{^)/^ for 
at Tb = 0.1. From the scaling form 
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S{ijj) — ^F{uj/Tb), the results for all Tb are recov- 
ered. We observed with very good precision that the 
"quasi-particle peak" j|] disappears at g = ^ (instead of 
the sometimes conjectured g = but we have no an- 
alytical proof of this result. Physically this means that 
the crossover between the oscillating behaviour and the 
damped one is at g = ^ at zero temperature. 

It is easy to check from the form factor expansion that 
for all g the following relation holds : 



lim ^LJ^ 



(20) 



which leads to the large time asymptotics C{t) « 

\/7r 1 

eg (7rTi3t)2 ■ 

At small frequency, Siuj) expands as a power series 
in [uj/TbY"' for any g. At large frequency, it expands 
as a double series in [TB/t^Y and {TbI^Y^"^^ ■ As a 
consequence, C(i) oc t^~'^^ at small times. 

In conclusion, we have shown that the description of 
quantum impurity problems based on (i) massless in- 
tegrable boundary field theories, (ii) boundary states, 
and (iii) form-factors, allows an efficient computation of 
the time dependent properties. So far, these had been 
completely inaccessible except by Monte Carlo simula- 
tions, and by various approximations which were recently 
proven unreliable |l^ . Our goal here was to present the 
technique and some results. Details will be provided in 
a more extensive paper. We also expect to be able to ex- 
tend the method to the case where there is a voltage (a 
bias) and to the finite temperature case. This will involve 
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a combination of both the form-factors and the TBA 
approach. Also, the careful reader might wonder why 
the computation of AC properties requires form-factors, 
while DC properties have been successfully computed so 
far by analogy with the free case. The reason is that DC 
properties see only the charge Q, ie the x integral of the 
current, and Q can be shown to act purely diagonally on 
multiparticle states. See for more details. 
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FIG. 1. Frequency dependent conductivity at T = 0, 
9 = 1/3. 

FIG. 2. Response function for different values of g. 
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